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ABSTRACT 

We present a multifrequency approach which optimizes the constraints on 
cosmological parameters with respect to extragalactic sources and secondary 
anisotropies contamination on small scales. We model with a minimal num- 
ber of parameters the expected dominant contaminations in intensity, such 
as unresolved point sources and the thermal Sunyaev-Zeldovich effect. The 
model for unresolved point sources, either Poisson distributed or clustered, 
uses data from Planck early results. The overall amplitude of these contribu- 
tions are included in a Markov Chain Monte Carlo analysis for the estimate of 
cosmological parameters. We show that our method is robust: as long as the 
main contaminants are taken into account the constraints on the cosmological 
parameters are unbiased regardless of the realistic uncertainties on the con- 
taminants. We show also that the two parameters modelling unresolved points 
sources are not prior dominated. 

Key words: Cosmology: cosmic microwave background - Physical data and 
processes: cosmological parameters. 



1 INTRODUCTION 

The CMB anisotropy data constitute a fundamental 
tool to test the standard cosmological model and its 
extensions. In particular small scale CMB anisotropies 
have a great importance in cosmology and are one 
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of the current frontiers in CMB observations. Many 
experiments have been dedicated to the observa- 
tions of small scale anisot r opies. The Planck satellite 
(|The Planck Collaboration! 120051 ) will cover all the 
scales up to £ ~ 2500 but there are also ground based 
experiments observing smaller regions of the sky but 
with higher angular resolutio ns, such as: ACBAR 



(|Reichardt et al l 120091). CBI dRedhead et al. | 12004), 
QUa D (IQUAD Collaboratiorj|2009l ), SPT (|Lueker et all 
l2010h . ACT jDas et al.ll201lh . At these scales, the Silk 
damping suppresses the primordial CMB contribu- 
tion with respect to extragalactic contamination and 
secondary anisotropies. These contaminations affect 
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the constraints on cosmological parameters derived 
from small scale data introducing biases and therefore 
it is necessary to properly account for them in the 
data analysis. We present an approach which has been 
developed to optimize, with respect to the contamina- 
tion by point source contributions, the constraints on 
cosmological parameters which could be obtained from 
high resolution multi-frequency data. 
The problem of extragalactic sources and secondary 
anisotropics contamination of the angular power 
spectrum and consequentely of the cosmological 
parameters has already been addressed in s evera l 
recent article s jDouspis. Aghanim fc Langerl 20061; 
Taburet et al.l 120081 ; iTaburet. Aghanim fc Douspij 



20091; lEfstathiou fc Grattonl 120091: iMillea et alj l201ll ; | Sunyaev fc Zeldovich 



MonteCarlo (MCMC) analysis for the cosmological and 
astrophysical parameters. In section 5 we present the 
discussion and conclusions. 



2 ASTROPHYSICAL CONTAMINATIONS 

One of the major contributions of astrophysical con- 
tamination on small scales, and in the regions of 
the sky clean from the Galactic emission, is given 
by discrete sources and in particular galaxy clusters 
and the extra-galactic point sources point sources. 
The former contrib utes through the Sunyaev Zel- 
dovich (SZ) effect jSunvaev fc Zeldovich I Il970.1980l ; 



19721 ) whereas the latter contri- 



Addison et alj l201ll ; lEfstathiou fc Migliacciol 12011 
together with the specific tools developed by the 
data analysis groups of the most recent small sca le 
experiments l|Keisler et al.ll201ll ; [Dunklev et alj|201ll ). 

The approach we present in this article consists 
in modelling the astrophysical contributions to the an- 
gular power spectrum with minimal parametrizations, 
where minimal relates to the number of parameters 
necessary to describe the signals, and considering the 
frequency dependences of the different contributions. 
This is the most impor tant difference with respect 
to previous approaches (lEfstathiou fc Grattonl 20091; 
Millea et all l201ll; I Addison et alj l201ll ; iKeisler et all 
201ll ; IDunklev et al.ll2oTli r These parametrizations are 
included as additional contributions to the primary 
CMB in the Markov Chain MonteCarlo analysis to 
constrain cosmological parameters. The parameters 
characterizing the astrophysical contributions are 
included in the analysis together with the cosmo- 
logical ones. As described in the following, we focus 
on the frequencies: 70,100,143,217,353 GHz. Wc 
use the available data fr om P lanck early results 
(|The Planck Collaboration! l2011al lbh to derive the 
parametrizations to describe the astrophysical signals. 
When no data are available, we rely on predictions 
from theoretical models and empirical simulations. 

In the following, we use a fiducial cosmologica l 
model with the parameters from IKomatsu et al.l (120111 ) 
which are reported in the third column of Table 4. 
The polarization from point sources and SZ effect is 
negligible we thus do not take it into account. We also 
do not account for residual signal, in intensity nor in 
polarization, from our Galaxy. 

The article is organized as follows. In Section 2 
we describe the astrophysical signals and we derive the 
parametrizations of their contributions to the angu- 
lar power spectrum. In Sect 3 we describe the analy- 
sis method and the frequency channel combinations. In 
Section 4 we present the results of the Markov Chain 



bution depends on the frequency. We have considered 
only the contribution of the thermal SZ effect since the 
kinetic SZ effect has a smaller amplitude. For the point- 
source contamination we considered different contribu- 
tions depending on the frequency. We aim in this section 
at parametrising the different contributions by reducing 
to their minimum the number of parameters necessary 
to characterize the signals, and by including the fre- 
quency dependence in these parameters. 

2.1 Thermal Sunyaev Zeldovich effect 

The thermal SZ effect (TSZ) is a local spectral distor- 
tion of the CMB given by the interactions of CMB pho- 
tons with the electrons of the hot gas in galaxy clusters 
( bunvaev fc Zeldovichl 11970.1 980 ) encountered during 
the propagation from the last scattering surface. 

The TSZ contribution to the CMB signal strongly 
depends on the cosmological model and the cluster dis- 
tribution through the dark matter halo mass function. It 
also depends on the cluster properties through the pro- 
jected gas pressure profile. The TSZ angular power spec- 
trum can be computed us i ng these two ingredients (e.g. 
IKomatsu fc Seliakl (|2002l); I Aghanim. Maiumdar fc Silkl 
(|2008l )). It was shown in lTaburet et alj (|2008l ) that the 
SZ contribution to the signal can be analysed in terms 
of a residual contribution to the total signal when the 
detected clusters are masked out from the maps, or 
in the contrary the TSZ signal can be analysed to- 
gether with the CMB one without taking out the de- 
tected clusters providing the total SZ power spectrum 
is pro perly modelled, see lTaburet. Aghanim fc Douspia 
(120091 ). The latter approach is the o ne that is used 
in th e CMB experiments lik e WMAP (IKomatsu et alj 
l201ll ) ACT ilDas et alj|201ll; IDunklev et al.ll201ll) and 
SPT jLueker et alj l2010l: IKeisler et al.ll201li) As dis- 
cussed in lTaburet. Aghanim fc Douspisl ( 20091 ) . we can 



consider three possible parametrizations for the TSZ 
power spectrum. However, the authors showed that 
the parametrization which represents physically the 
best the TSZ signal and at the same time does not 
bias the cosmological parameters is the one derived in 
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iKomatsu fc Seliakl (2002). It considers a template for 
the spectral shape, Ct, and a normalisation which de- 
pends on the parameters as and Q^h 2 . According to 
these results, we chose to use: 



Cf (y) = A sz g(u)aln 2 b h 2 C e , 



(2.1) 



where Asz accounts for the uncertainty in the normal- 
isation. The frequency dependence of the signal for the 
TSZ is given by the analytical function g(v jj. In the fol- 
lowing, we used a fi ducial TSZ power-s pectrum shape 
computed using the lTinker et all J200ST) mass function 
and a pressure profile provided in lArnaud et al.l (|2009h . 



2.2 Point source contribution 

An experiment like Planck with its moderate angular 
resolution is not optimized for point source detection. 
However, Planck, thanks to its full sky coverage and 
its 9 frequencies, is detecting a large number of sources 
with the d etection thresholds reporte d in Ta ble [T] and 
taken from iThe Planck Collaboration! i|2011rj ). For the 
CMB analysis, the contribution of detected sources 
is removed by masking them in the maps. However, 
the contribution of unresolved point sources remains 
and needs to be properly dealt with since it acts like 
an unavoidable noise that impacts the angular power 
spectrum amplitude and the cosmological information 
contained in it. 

In the frequency range of interest for our study the 
main contributions from extragalactic point sources are 
those of the radio galaxies which dominate the low fre- 
quency channels and the infrared ones which dominate 
the higher frequencies. Together with these two main 
populations there is a contribution from anomalous 
objects like Gigahertz Peaked Spectrum, Advection 
Dominated Accretion Flows/ Adiabatic Inflow Outflow 
Solution and radio afterglows of Gamma Ray Burs t 
(|Toffolatti et al.lll998l . 120051 ; |Pe Zotti et ai1l2005l . 120061 ) . 

Radio galaxies, typically BL LACs and Flat Spec- 
trum Radio Quasar s present a f l at sp ectrum 5V oc v a 
with a > —0.5 l|Danese et al. 1987). There is also 
a contribution of steep spectra radio sources which 
typically have fluxes S v oc i> a with a ~ —0.7 — (—1). 
The anomalous objects have only a minor contribution 
to the total of ra dio galaxy power at the arcmin scale 
(see for example iToffolatti et al.l l|2005l ) for detailed 
discussion). Infrared galaxy emission presents a steep 
spectrum S v oc v 3 ' 5 and is due to the UV light 
from star formation reprocessed by dust in galaxies. 



1 In the non-relativistic approximation, it is given by: g(i>) 
~ (*fS=T " 4 ) With X = 56.78GHz - 



Infrared galaxies have been lengthly studied and ob- 
served in the recent past especially by CM B-dedicated 



experiments like ACT 
ences therein), SPT 



T (iDas et all J201ll) 
dLueker et alj l20ld ) 
ationlbOllbT) but alsc 



The Planck Collaboration! 2011b) but also by S pitzer 



and refer- 
and Planck 



Ashbv et al.l 120091^ . Blast (iMarsden et all \200$ ). and 
Herschel (|Amblard et alj|201 ll) . 



In our range of frequencies and for our considered 
thresholds, the radio source contribution to the angu- 
lar power spectrum can be described only by a Poisso- 
nian term. The situation is quite different for infrared 
galaxies. The latter have in fact a strong clustering con- 
tribution to the angular power specturm and a care- 
ful modeling and p arametrisation of their clustering 
term are necessa ry (|The Planck Collaboration] l2011bl : 
IPenin et al.ll201lh . We thus consider the following con- 
tributions from point sources: 

70 GHz -> Poisson(Radio) 

100 GHz -> Poisson(Radio) + IR Clustering 

143 GHz -> Poisson(Radio + IR) + IR Clustering 

217 GHz -> Poisson(Radio + IR) + IR Clustering 

353 GHz -> Poisson(Radio + IR) + IR Clustering . 

In the following, we describe the different contri- 
butions from point sources to the CMB anisotropy, we 
detail the proposed parametrisations and then provide 
the associated angular power spectrum. 

2.3 Poissonian contribution 

The Poisson contribution of point sources is given by 
their random distribution in the sky. The coefficients 
of spherical harmonics of a Poissonian distribution of a 
jopulation of sources in the s ky with flux density S are 



population 01 sourc es in the s k 
(|Tegmark fc Efstathiou1 ll996): 



(ai m ) 



V^nnS for £ = 
for f/0,Vm 



where n = N/Ati is the mean number of sources per 
steradian with flux density S. The angular power spec- 
trum is given by Ce = (\aim\ 2 ) — |(ei£ m )| 2 = nS 2 , which 
can be generalized to Ci = E^Sf if we consider sources 
with different flux densities. Extrapolating the contin- 
uum limit in flux up to the detection threshold, S max , 
the angular power spectrum for a Poissonian distribu- 
tion of sources in sky is given by a flat Cf. 



.^-tPois 



dJV(g) s 



(2.2) 



where di ^^ are the source number counts. The flat 

dS 

spectral shape of the Poissonian contribution makes 
it easy to parametrize: It requires a simple amplitude 
term. One possible choice of parametrization is to use a 
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different free parameter for the amplitude for each fre- 
quency as in lEfstathiou fc Grattonl (|2009h . This would 
require in our case for the radio sources alone 5 param- 
eters for the Poisson term. Since we aim at reducing the 
number of parameters to the minimum, we construct a 
single amplitude parameter for the total Poisson con- 
tribution and accounting for the frequency dependence. 
This frequency dependence is complicated by the de- 
pendence on the detection threshold S max and thus no 
analytical expression can be constructed. We therefore 
reconstruct from the data and from the galaxy evolution 
models the necessary information to derive the single 
parameter for the Poisson term. 

Both radio and infrared galaxies contribute with a 
Poissonian term. We treated the former by fitting the 
number counts and then integrating on the flux den- 
sity whereas for the latter we used the available results 
from (Bethermin et al.ll201ll ; IThe Planck Collaboration! 
12011b ). The two contributions are then summed in a 
single Poissonian term. The expected Poissonian con- 
tribution for each frequency can be computed using the 
number counts given a detection threshold. We used 
the predi cted number counts for rad iogalaxies from the 
model of (|De Zotti et al.ll2005l , [2fJ06l ). It is a good repre- 
sentation of low flux part of the observed number counts 
but the comparison with the Planck Early Release Com- 
pact Source Catalogue (ERCSC) showed a bias at high 
flux densities with t he De Zotti et al. model over- 
predicting the co unts l|The Planck Collaborationll2011al : 
iTucci et al.ll201ll ). Since the ERCSC number counts arc 
limited to bright flux densities, we constructed com- 
posite number counts combining the De Zotti model 
predicted counts, for lower flux densities roughly below 
0.2 — 0.3 Jy (thin dots in Figs. [T] and [2]), whereas for flux 
densities greater than 1 Jy we used the actual Planck 
measur ements reported in IThe Planck Collaboration! 
(^OllaT ) (thick dots in Figs. [Q and [2|. 

At 353 GHz, the radio sources are subdominant 
with respect to infrared sources. In order to compute 
the associated Poisson term for the 353 GHz channel 
we therefore used the predictions of the De Zotti et al. 
model that we rescaled to account for its likely overpre- 
diction, similar to that at 143 and 217 GHz. In practice, 
we have extrapolated the correction factors for the ra- 
diogalaxy Poissonian contribution to the angular power 
spectrum of the 143 and 217 GHz channels, 0.93 and 
0.92 respectively, and found a correction factor of 0.91 
for the Poissonian angular power spectrum of the 353 
GHz channel. 

The obtained number countfl are represented for 
each frequency by the weighted sum of three differ- 
ent contributions: one at low i 7 "^ ", one at intermediate 
F™ ed , and one at high F^ lgh flux densities. The number 



2 The differential number counts are normalized to the Eu- 
clidean ones ( di ^i^ <x S~ 2 ' 5 ). 



counts are given by 
2 . 5 dN(S) tot 



S 



dS 



Ei = Low, Med, High wl Fl 



(2.3) 



Each term i^°w,Med,High = S -2.5diV|S) writeg ag 
F L OW ,Med,m g h = For each term and each 

frequency, the exponents a and /3 and the coefficients 
A, B and C, together with the weights of the sum 

w ^ow,Med,High &r<} ^ ad j usted by compar i n g them 

to the composite number counts. Table 5 summarises 
the set of parameters used to construct the composite 
source counts and the fit to the hybrid model/data 
counts is shown as the thin line in Figs. [1] and [2] We 
note the very good agreement between the analytical 
fit and the theoretical/measured points. 

In Figs. [3] and [4] we compare o ur hybrid number 
counts to those of ITucci et al.l l|201ll ) based also on the 
ERCSC. Our fit overpredicts the number counts in the 
range 0.01 — 1 Jy with respect to the Tucci et al. model. 
This overprediction is due to the fact that the flux den- 
sity at we match the predictions and the data is 0.01 Jy 
whereas Tucci et al. chose a lower value. 

Using the obtained composite counts we compute 
the Poissonnian power spectra at the considered fre- 
quencies. They are given by following expressions: 



'max, i> 



x 2F\ [cu, „, 6j, u,Ci,v, di, ^£S a x, v MM) 



where 2F1 are the Hypergeometric functions of sec- 
ond type ijAbramowitz fc Stegun|["l965r i. the sum on i is 
over the three contributions Low, Med and High, and 
the coefficients and exponents depend on the frequency 
and are reported in Table 6. The values of the Poissonian 
contribution to the angular power spectrum that we ob- 
tained from our fitted number counts are reported in the 
third column of Table [T] For comparison in the fourth 
column we report the values computed from the origi- 
nal De Zotti-model fits. As expected from the number 
counts comparison, we have a higher Poissonian con- 
tribution than the one obtained with the Tucci et al. 
model, namely 30% higher at 70 and 100 GHz, 40 and 
60% at 143 and 217 GHz respectively. At 353 GHz we 
predict an amplitude twice larger. 

Together with the term coming from the ra- 
dio sources it is necessary to consider the Poisso- 
nian contribution of infrared galaxies which is rel- 
evant only for the higher frequency channels the 
143, 217 and 353 GHZ. The values of the Poisso- 
nian infrared contr ibutions at these frequ encies are 
computed from the iBethermin et al.l {2011) model in 
IThe Planck Collaboration! (|2011bT > and reported in the 
fifth column of Table [T] We summed the two Poisson 
terms for radio and infrared in the frequencies in a sin- 
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Figure 1. Fits of the hybrid model/data number counts for 
each frequency. The solid line is the empirical fit whereas 
the dots represent the De Zotti et al. model number counts 
(thin light dots) and Planck data (dark thick dots). The three 
curves under the fit are the three contributions (Low, Med 
and High, respectively dotted, dashed and dot-dahsed lines) 
which sum in the total fit. From top to bottom: 70 GHz, 100 
GHz, 143 GHz. 



gle Poissonian term: 



C 



Apse™ 



PS,Radio+IR 



(2.5) 



where C t „ is the Poissonian contribution at the different 
frequencies whereas Aps which incorporates deviations 
of real data from this model remains the same for all 
the frequency considered. 



o.oi - 

0.002 V, ' 



1.x 10~'l.x KT'O.OOOOl 0.0001 0.001 0.01 0.1 1 10 

S(Jy) 

353 GHz 




I.xl0"'l.xl0-"0.00001 0.0001 0.001 0.01 0.1 1 10 

S (Jy) 




dO-'O.OOOOl 0.0001 0.001 
S(Jy) 

Figure 2. Fits of the hybrid model/data number counts for 
each frequency. The solid line is the empirical fit whereas 
the dots represent the De Zotti et al. model number counts 
(thin light dots) and Planck data (dark thick dots). The three 
curves under the fit are the three contributions (Low, Med 
and High, respectively dotted, dashed and dot-dahsed lines) 
which sum in the total fit. From top to bottom: 217 GHz, 353 
GHz. Bottom is the comparison of the fits for the different 
frequencies: 70 GHz is purple, 100 GHz is cyan, 143 GHz is 
green, 217 GHz is yellow, 353 GHz is orange. 



2.4 Contribution from clustering term 

The Cosmic Infrared Background (CIB) is the inte- 
grated IR emission over all re dshifts of unresolved in- 
frare d star forming galaxies (|Lagache. Puget fc Dole! 
2005). This contribution is mainly that of star- 
bust galaxies, Luminous InfraRed Galaxies (LIRGs, 
IO^Lq < Lir < 10 12 L Q ) and Ultra Luminous In- 
fraRed Galaxies (ULIRG s, Lir > 10 12 Lp)), (with small 
contribution from AGN ) (Lagachc, Dole & Pugrt j [20031; 
lLagache. Puget fc Dole! 120051 : Fernandez-Conde et alj 
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70GHz 




0.00001 0.0001 0.001 0.01 0.1 1 
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100GHz 
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s 



143GHz 




0.00001 0.0001 0.001 0.01 0.1 1 



s 

Figure 3. Comparison of the fits of the hybrid model/data 
number counts with the Tucci et al. model for each frequency. 
The solid line is the empirical fit whereas the cyan dots rep- 
resent the De Zotti et al. model number counts (thin light 
dots) and Planck data (dark thick dots). The Tucci et al. 
model is represented by pink dots. From top to bottom: 70 
GHz, fOO GHz, 143 GHz. 



217GHz 




0.00001 0.0001 0.001 0.01 0.1 1 
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353GHz 
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Figure 4. Comparison of the fits of the hybrid model/data 
number counts with the Tucci et al. model for each frequency. 
The solid line is the empirical fit whereas the cyan dots rep- 
resent the De Zotti et al. model number counts (thin light 
dots) and Planck data (dark thick dots). The Tucci et al. 
model is represented by pink dots. From top to bottom: 217 
GHz, 353 GHz. 



Channels 


5Vnax 


.^Pois, Radio 
^Hybrid 


£tPois, Radio 
Z Original 


<^<Pois 
U £ IR 




(Jy) 


(MK 2 ) 


(MK 2 ) 


(^K 2 ) 


70 


0.57 


1.13 x 10~ 3 


1.08 x 10" 3 




100 


0.41 


2.05 x 10" 4 


2.25 x 10~ 4 




143 


0.25 


4.52 x 10~ 5 


4.86 x 10 -5 


1.2 ± 0.2 x 10~ 5 


217 


0.16 


1.60 x 10~ 5 


1.74 x 10 -5 


6± 1 x 10~ 5 



353 0.33 9.40 x 10~ 5 1.03 x 10~ 4 1.9 ± 0.3 x 10~ 3 



2008). One of the complexity related to the CIB is 
that at different frequencies it accounts for contribu- 
tions of IR sources from diffe rent redshift ranges (e.g. 
iFernandez-Conde et all i|2008l )). 

The c lu sterin g power spectrum is g iv en b y 
Knox et all d200 it) : IFernandez-Conde et all i|200St ); 
Penin et all |201ll ): 

Cfust = [ dz ?M^ x 
J d A dz 

jl(z)P g (k)\ k= t /dA (z) , (2.6) 
where r is the proper distance, cLa is the comoving angu- 



Table 1. Poissonian contributions to the angular power spec- 
trum for each frequency. In the first column are reported 
the observational frequencies and in the second the detection 
threshold fluxes. In the third column are reported the values 
of the Poissonian angular power spectrum for radio sources 
obtained using the composite number count fits. In the fourth 
column are reported the values of the Poissonian angular 
power spectrum computed using the original De Zotti et al. 
model number counts instead of the composite ones. In the 
last column are reported the value s of the Poissonian term for 
infrare d sources as predicted in l|The Planck Collaboration! 
l2011bh . 
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lar diameter distance, k = £/<1a derives from the Limber 
approximation, the Jr^§ C[2 (- z ) term takes into account 

A 

all geometrical effects and depends on the cosmological 
model adopted, j\(z) is the mean galaxy emissivity per 
unit of comoving volume and P g is the galaxy three- 
dimesional power spectrum . The emissivity of infrared 
galaxies can be written as l|Penin et al.ll201lh : 



•<"-(•£) 



S(Lj 



dN 



dzd(lnLi 



-d(lnLiR) (2.7) 



wnere dzd(inL IR ) * s ^ ne l ummos ity function and S is 
the flux. The clustering term as shown by Eq. 12.61 de- 
pends on the cosmological model through the geomet- 
rical factor, and through the galaxy power spectrum. 
Another complication is given by the indirect depen- 
dence of the emissitivity on the as sumptions on the 
galaxy Spectral Energy Distr ibution |Knox et al.ll200ll ; 
iFernandez-Conde et alj[2008h . Finally, the bias, relating 
the the emissivity fluctuations to the dark matter den- 
sity fluctuations, also affects the clustering term and is 
poorly known. In view of the the theoretical complexity 
of the clustering term, we rather used the observational 
clustering term extr acted directly from Planck maps a t 
217 and 353 GHz bv lThe Planck Collaboration! (|2011bh . 
In order to minimally paramatrize this contribution, we 
have fitted the Planck best fit models. Since the Planck 
data are not available for the channels below 217 GHz, 
the values for the clustering term for 100 and 143 GHz 
channels were derived using the halo occupation distri- 
bution b est-fit parameters of the 21 7 GHz c hannels (Ta- 
ble 7 of iThe Planck Collaboration! <|2011bh ), combined 
with th e emissivities at 143 and 100 GHz computed us- 
ing the iBethermin et al.l (|2011T ) model. 

The fit to the clustering term is given for all the 
frequencies by the following expression: 



£{£+!) 



c; 



clust 

_e 

2tt 



(v) oc 



logW 



(100 + £)tM 



,(2.8) 



where <f)(£,v) is a fifth order polynomial term Ao/£ 2 + 
Axil + A 2 + A 3 £ + Ad 2 + A 5 £ s + A 6 £ 4 accounting for 
the steepening at high £, whereas the low £ shape is 
dominated by a complex logarithmic form. The coeffi- 
cients Ao-e depend on the frequency. The first step to 
obtain a fit to the clustered term was derive a spectral 
shape, namely the combination of a polynomial compo- 
nent with the logarithm form. Then we fitted the form 
given in Eg. 12.81 for each frequency first fitting the loga- 
rithmic function and then the polynomial function. Fi- 
nally, we have derived the frequency dependence of each 
coefficient Ao-e and each exponent by fitting to the dif- 
ferent frequency. The resulting fit to the power spectrum 
is given by: 

log(<0 A » M 
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Figure 5. Comparison between the best fit models from data 
(for 217 and 353) and extrapolations (for 100 and 143) of the 
infrared source clustering (dashed) and the fits (solid line) 
for the various frequency channels. Colors: Cyan 100 GHz, 
Green 143 GHz, Yellow 217 GHz, Red 353 GHz. 



where the first sum, on i, is on the coefficients of the 
polynomial and the second, in j is on the polynomial de- 
pendence on the frequency. In Table 7 are reported the 
coefficients a\. Such complex formula insures an accu- 
racy of the fit at < 1% level in all £ range, and at < 0.1% 
in about 90% of the considered £ range. 

Finally, provided the expression Eq. 12.91 giving the 
frequency-dependent CIB clustering power spectrum, 
the parametrization of the clustering term can be simply 
written as: 



Cf us \v) = ylcLCf ust H 



(2.10) 



where the amplitude parameter Acl accounts for devi- 
ations of the real signal from our parametrisation. 



3 METHODOLOGY 

In the previous section we have detailed how we con- 
struct a frequency dependent parametrization for the 
point sources and TSZ contributions to the angular 
power spectrum. In Fig. [6] we display each contribution 
together with the primary CMB spectrum for the five 
frequencies of interest. As anticipated, we note that the 
higher frequency channels are dominated by the cluster- 
ing term contribution, the 353 GHz channel being the 
most contaminated. The lower frequencies are in turn 
dominated by the Poissonian term contribution. We also 
note that at small scales the TSZ contribution remains 
subdominant with respect to that of the point sources. 
As seen in Fig. [6] displaying the comparison between 
the sum of all contributions with the primary CMB, 
the 353 GHz channel is the most contaminated due to 
the clustering contribution, whereas the cleanest chan- 
nel is the 143 GHz due to the low contribution of both 
clustering and Poissonian terms. We show in Table [2] 
the relative importance of the two main contributions 
of clustering and Poissonian expresssed in terms of the 
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Channels 


^Clustering 


^.Pois, Radio 


CT Pois,IR 






Jy/s/sf 


Jy/Vsr 


100 


193.906 


789.309 




143 


427.572 


786.607 


405.302 


217 


2890.06 


873.219 


1690.98 


353 


14927.1 


1312.22 


5899.56 



Table 2. Rms of the different contributions compared. In 
the first column are reported the interested frequencies. In 
the second column are reported the values of the clustering 
rms, in the third the Poissonian rms for radio sources and in 
the fourth the rms values for the Infrared Poissonian. 



root mean square: 

a = v /E ; (2Z + l)6fC ; /(4 7 r). (3.11) 

We note again how the high frequencies have a dominant 
clustering contribution whereas for the 100 and 143 GHz 
the dominant contribution is the Poissonian term. 

Since we are interested in small scale anisotropics 
it is possible to avoid the component separation process 
to remove diffuse emission from our galaxy and use in- 
stead single frequency angular power spectra derived by 
masking the Galaxy and the resolved point sources. The 
Galaxy may contribute even after masking with high 
latitude emissions and residuals which may remain de- 
pending on the aggressiveness of the mask used. In this 
work, we neglect the possible residuals, after masking, 
from Galactic contamination. 

We combine the angular power spectra of the mi- 
crowave sky at the considered frequency channels in a 
single power spectrum that we use to estimate the cos- 
mological parameters. An optimal method to combine 
channels is the inverse noise variance weighting scheme 
(jEfstathiou fc Grattonl2009l '). This weighting scheme in- 
deed gives the larger weights to the channels with higher 
resolutions and lower noise levels. This scheme is par- 
ticularly adapted to the case of high multipole analysis 
where the noise and the beam smearing start to domi- 
nate. We therefore use a weighted linear sum of single 
frequency power spectra: 

c\ ot = V iWi (l)Cl 

where the summation in i is over the frequencies. 

For simplicity, the noise is modelled as an isotropic 
Gaussian white noise with variance a% on the total in- 
tegration time. We define the beam function as: 

y2 _ e -«(£+l)(0.425FWHM/60 7r/180) 2 . ^ ^) 

where FWHM is the full width half maximum of 
the channel-beams in arcminutes. We illustrate our 
method on a Planck-like c ase taking the beams 
and noise characterist i cs fro m iPlanck HFI Core Team! 
(|201lMMennella et all (|201ll ). Moreover, the noise has 



Instrument 


LFI 




HFI 




Center frequency GHz 


70 


100 


143 


217 


353 


Mean FWHM (arcmin) 


13 


9.88 


7.18 


4.71 


4.5 


AT per pixel (7) 


24.2525 


8.175 


5.995 


13.08 


54.5 


AT per pixel (Q kU) 


34.6075 


13.08 


11.1725 


24.525 


103.55 



Table 3. Planck channel performance charac teristics 
jMennella et alj|201ll : TPlanck HFI Core Teamll201lf) for the 
nominal mission. 



been divided by 2 to account for the approved exten- 
sion of the mission. The numbers are summarised in 
Table El 

We can thus define the noise function as: 



(<7jv) f2pi x 
b? 



(3.13) 



where O p i x is the area corresponding to the pixel. Note 
that the total noise contribution to the angular power 
spectrum considers also the co ntribution of the observed 
sky fraction AQ/Cl = N e y/2/(f aky (2£ + 1)). Assum- 
ing that the sky fraction in the different channels is 
the same, the contribution of this term cancels in the 
weights. With these definitions the weights for the single 
frequency spectrum is given by: 



Wi(£) 



v. 1 

J W 



(3.14) 



where i stands for the specific channel whereas j runs 
over all channels. 

We show the weights obtained using the inverse 
noise variance scheme in Fig. [7] We note how the dom- 
inant contribution among the various combinations is 
given by the 143GHz and 217 GHz channels which have 
higher resolution and lower noise level, whereas channels 
with lower angular resolution have very little weights 
like the 70 GHz whose contribution is strongly subdom- 
inant within this approach. 



4 ESTIMATING THE COSMOLOGICAL 
AND EXTRAGALACTIC PARAMETERS 

We now show the results of our analysis where we have 
applied, to the mock data, our approach of frequency- 
dependent minimal parametrisation of the extragalactic 
sources. We first show the capability of our approach to 
perfectly recover the input cosmological parameters in 
presence of extragalactic contamination and exhibit the 
importance of the wide frequency range for the charac- 
terization of the extragalactic signals by comparing the 
results for different frequency combinations. Then we 
show the impact of the astrophysical contributions to 
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Figure 6. Comparison of the extragalactic signals with pri- 
mary CMB (black solid line) for the different frequencies. 
The line color represent the frequency: 70 GHz (purple), 100 
GHz (cyan), 143 GHz (green), 217 GHz (yellow), 353 GHz 
(red). In the upper panel we show the comparison for each 
frequency of the three extragalactic signals: dotted is the 
clustering term, dashed is the Poissonian term and solid is 
TSZ. In the lower panel we show the comparison of the total 
astrophysical contribution with the CMB for each frequency. 
The solid lines are the sum of the TSZ, Poissonian and clus- 
tering contribution and follow the color scheme above. 
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Figure 7. Weights of the different frequencies with the in- 
verse noise variance scheme: purple 70 GHz, cyan 100 GHz, 
green 143 GHz, yellow 217 GHz, red 353 GHz. 



the power spectrum and to the cosmological parameters 
and their errors. We finally test the robustness of our 
approach by investigating the dependence of the results 
on the astrophysical models we used. We will show 
also the results in the extended parameter space in 
particular we considered the case of a ACDM with 
running of the spectral index and tensor modes. 

We developed an extension of the pub licly avail- 
able CosmoMC code |Lewis fc Bridle! l2002t ) which, in 
addition to estimating the cosmological parameters, es- 
timates at the same time the contributions of a Poisso- 
nian term, a clustering term and a TSZ term. This is 
done by deriving a minimal number of amplitude terms, 
namely Asz , Aps , Acl , which represent the departures 
from the input point source and TSZ models. 

In practice, we have modifie d the CAMB code 
(|Lewis. Challinor fc Lasenbvl I2OO0I ) . version May 2010, 
to compute for each given frequency channel the sum 
of point-source and TSZ spectra as provided by our 
parametrization (Sect. 2) and the standard CMB angu- 
lar power spectrum. The resulting total power spectra 
at the different frequencies are then combined in a single 
effective power spectrum in the CosmoMC code. 

We first explore the frequency combination and its 
effect on the recovery of the cosmological and the astro- 
physical parameters. The latter are represented by the 
amplitude^] Asz, Aps and Acl- We have considered the 
results using five combinations of frequencies: 143 and 
217 GHz, 100 to 217 GHz, 100 to 353 GHz, 70 to 217 
GHz, and finally all the frequencies under consideration 
70 to 353 GHz, with a sky coverage of / s k y = 0.85. In 
order to establish which frequency combination recovers 
the cosmological and astrophysical parameters best, we 
have performed for each combination an MCMC anal- 
ysis of mock data using a ACDM-based CMB plus the 
astrophysical contributions. For the latter we used the 
"nominal" contributions, which means all the amplitude 
terms were set to one (Asz = Aps = Acl = 1). All the 
channel combinations perfectly recover the cosmological 
parameters. The major differences are in the recovery 
of the astrophysical parameters, Asz, Aps and Acl, see 
Fig. [8] More specifically, the simplest combination of 
the 143 and 217 GHz channels recovers the Poissonian 
term with the smallest errors but it gives large errors 
for the clustering term. By contrast, the combination of 
all considered channels, 70 to 353 GHz, gives the largest 
errors for the Poissonian term but significantly reduces 
the errors on the clustering term. The clustering term 

3 Note that our treatment of the TSZ contamination samples 
on all the physical factors multiplying the template Ci in Eq. 
12.11 i.e. Asz , f 8 , Of, , h. Therefore our approach differs from 
that adopted by the WMAP team which samples only on the 
overall amplitude multip lying the Komatsu-Seljak template 
l lKomatsu fc SeliaJj|200a) . 
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Figure 8. Comparison of the results for the two main astro- 
physical parameters of the analysis obtained with different 
frequencies, always combined with the inverse noise variance 
scheme: black (solid) is 70-353 GHz, red (short dashed) is 
143-217 GHz, green (long dashed) is 100-217 GHz, blue (dou- 
ble dot dashed) is 70-217 GHz and the pink (dotted) line is 
the 100-353 GHz. 



Parameter Best Fit value I Best fit value II Input value 



n b h 2 


0.02251 ± 0.0001 


0.02254 ± 0.0001 


0.0225 


n c h 2 


0.1121 ± 0.001 


0.1122 ±0.001 


0.112 


T 


0.0885 ± 0.004 


0.0874 ± 0.004 


0.088 


n 3 


0.967 ± 0.003 


0.966 ± 0.004 


0.967 


log[A s 10 10 ] 


3.19 ± 0.007 


3.19 ±0.007 


3.19 


H 


69.96 ±0.5 


69.94 ± 0.54 


70 


A PS 




0.99 ± 0.17 


1 


A CL 




1.01 ± 0.22 


1 


A S z 






1 



Table 4. Comparison between the resulting best fit and the 
input values for the cosmological and extragalactic parame- 
ters. The fiducial model for the mock data is reported in the 
fourth column. In the second column are given the parame- 
ters recovered when taking into account the nominal contri- 
butions of the point sources and TSZ signals (note that the 
Asz remains unconstrained). In the third column the pa- 
rameters recovered when varying also the parameters of the 
contributions of the point sources and TSZ signals. The CMB 
plus extragalactic contribution power spectrum obtained us- 
ing the best fit values and the one generated using the fiducial 
model values do not differ by more than 0.2%. 



being the most uncertain with respect to the Poissonian 
one, we chose to use the combination of the five channels 
in order to improve the constraints on this major con- 
tribution to the microwave sky at small angular scales. 

4.1 Standard ACDM model 

We performed the analysis in the framework of a stan- 
dard ACDM cosmological model. The input mock Ct 
were produced with nominal values for the point source 
and TSZ contributions, namely Asz = Aps = Acl = 
1. Two cases are compared: Fixed values for the as- 
trophysical parameters (fixed to one), and free values 
for ylsz, ^4ps, and Acl- The results are displayed in 
Fig. and summarised in Table 4 (second column is 
the case with astrophysical parameters fixed to one, the 
third column is the case with astrophysical parameters 
variation), together with the input parameters (fourth 
column) . We show that the input cosmological param- 
eters are accurately recovered in both cases. When the 
astrophysical parameters Asz, Aps, and Acl together 
with the cosmological parameters are let free the er- 
rors bars are slightly enlarged. The astrophysical pa- 
rameters, Aps and Acl, which characterize the point 
sources (Poissonian and clustering terms) are recovered 
well but they exhibit a strong degeneracy. This degener- 
acy is expected, in fact the Poissonian and the cluster- 
ing terms have different shapes but they both contribute 
additively to increase the measured power at high mul- 
tipoles. The results are the same when we use the Tucci 
et al. model rather than out hybrid model for the radio 
contribution. The amplitude of TSZ spectrum in turn 
remains unconstrained due to its subdominance with 



respect to the other two contributions of point sources. 
We will thus not discuss it further. 

We now explore the robustness of the approach, in 
terms of parameter recovery, with respect to the mod- 
els we use for the point source contributions. We do not 
vary the TSZ signal as it is subdominant and cannot be 
recovered. For the effects of TSZ on the cosmological pa- 
ra meters in absence of p o int sources, we refer the reader 
to iTaburet et~ai] i|2008r ); iTaburet. Aghanim fc Douspisl 



(|2009l ). We generated mock data including the CMB, 
the nominal TSZ contribution and the point source 
contributions (Poisson and clustering terms). The later 
have been chosen to differ significantly from the nominal 
contributions. 

We first focus on the Poissonian term and multi- 
plied the nominal Poissonian term by two in the mock 
data for each frequency channel leaving untouched the 
other contributions. In Fig. HOI we present the obtained 
cosmological and astrophysical parameters. It is worth 
noting that even with a much larger Poissonian contri- 
bution the cosmological and the astrophysical parame- 
ters, Aps and Acl, are perfectly recovered. 
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Figure 9. Cosmological (and astrophysical) parameters with and without varying the astrophysical contribution parameters. 
The five frequency channels are combined with the inverse noise variance weighting scheme. The black (solid) curve represents 
the analysis which in addition to cosmological parameters varies also the astrophysical ones whereas the red (dashed) curve 
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Figure 10. Robustness test part 1. Cosmological and astro- 
physical parameters with different astrophysical signals with 
respect to the fiducial model. The black (solid) curve are 
the results of the analysis which considers a Poissonian term 
which is two times its fiducial value. The red (dashed) curve 
represent the analysis with a clustering term for the 100 and 
143 GHz doubled with respect to the fiducial value. 



If we now focus on the clustering term, many tests 
with respect to the nominal model can be envisaged 
to test the effect of the parametrisation on the pa- 
rameter recovery. We test the robustness of the re- 
sults, or in other words the sensitivity to the cluster- 
ing term, with respect to the frequency dependence 
of the clustering term by varying not only the ampli- 
tude but also the frequency dependence of the cluster- 
ing term. The clustering term at 353 GHz being very 
well constrained by the present Planck measurements 
of the CIB, we choose not to vary it. The CIB measure- 
men ts at 217 GHz by Planck, in turn, h ave larger errors 
(see lThe Planck Collaboration! l|2011bl )) and the values 
adopted in the parametrisation for the 100 and 143 GHz 
channels are empirical extrapolation based on the 217 
GHz measure. Therefore it is very likely that at these 
frequencies the clustering term may be significantly dif- 
ferent from what we used in the model. For this reason, 
we chose to vary these three frequencies. 

We first analyse a case where we multiply by two 
the clustering contribution of the 100 GHz and 143 GHz 
channels. In Fig. [10] the results show that this modifi- 



cation does not affect significantly the cosmological nor 
the astrophysical parameters. The reason is that at 100 
and 143 GHz the clustering term is not the dominant 
astrophysical component and therefore in the frequency 
combined spectrum the variation of the clustering term 
for these two frequencies does not modify the final re- 
sult. In the inverse noise variance weighting combination 
scheme the major contribution to the clustering term 
comes from the 217 GHz channel which has a strong 
clustering component and at the same time is the domi- 
nant channel at high multipoles in the frequency combi- 
nation. Since the 217 GHz measurements of the cluster- 
ing still allow a rather significant variation, we decided 
to test the effects of a variation of the 217 clustering 
contribution by varying it by a factor 2. The obtained 
contribution differs by several sigmas from the present 
constraints obtained by Planck. We see from Fig. llll that 
the cosmological parameters are perfectly recovered ei- 
ther when we multiply or when we divide the 217 GHz 
clustering contribution by two. Unsurprisingly, the re- 
covered mean values of Acl are multiplied/divided by a 
factor near two following to the varied contribution. Due 
to the degeneracy between the Poissonian and cluster- 
ing terms, the mean values of Aps are also shifted with 
respect to the input value. 

We now test the robustness of the cosmological pa- 
rameter recovery with respect to the frequency depen- 
dence that we proposed in our parametrisation of the 
clustering term. To do so, we multiply by two the clus- 
tering terms of the 100 and 143 GHz, whereas the clus- 
tering term at 217 GHz is multiplied and then divided 
by two. Here again the cosmological parameters are per- 
fectly recovered in both cases whereas the amplitudes of 
the Poisson and clustering terms vary. The results ob- 
tained for the mean values of Aps and Acl are the same 
as those resulting from varying only the 217 GHz con- 
tribution. This comes from the dominance of the 217 
GHz channel weight, in the channel combination, which 
in turn imply a dominance of the clustering term in the 
analysis. 

We also tested the case where the Poissonian contri- 
bution is completely different from our model. To do so, 
we analyse the data using our hybrid model whereas we 
construct the mock data from the Tucci et al. model. By 
doing so, the Poissonian term varies from frequency to 
frequency whereas the IR contribution is the same as in 
the previous tests. In Fig |12l we show that this different 
Poissonian contribution does not alter the cosmological 
parameters but the astrophysical ones are modified. In 
particular the Poissonian contribution moves along the 
degeneracy line and the clustering amplitude is shifted. 

4.2 Extended cosmological parameter space 

We have shown the results of the cosmological and astro- 
physical parameter recovery for a standard Lambda six- 
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Figure 11. Robustness test part 2. Cosmological and astro- 
physical parameters with different astrophysical signals with 
respect to the fiducial model. The black (solid) curve are the 
results of the analysis which considers the 217 GHz clustering 
term double its fiducial value. The red (dashed) curve repre- 
sent the analysis with a clustering term for the 217 GHz half 
its fiducial value. 



parameter cosmological model. We now present some re- 
sults obtained in the case of an extended cosmological 
parameter space. 

We considered together with the six standard pa- 
rameters both tensor perturbations and the running of 
the spectral index, varying the tensor-to-scalar ratio, r, 
and the running, n Iun , together with the six standard 
parameters and the three astrophysical ones. We have 
generated the mock data with both extra parameters, 
n run and r set to zero. We fix the tensor spectral index 
to the second order inflationary consistency condition 
n T = (2 - 4f - ns) ■ We included only the TE and 
EE mode in polarization we do not include the B-mode 
polarization power spectra in the exact form of the CMB 
likelihood we exploit, since the main limitation in that 
respect is given by the ability of removing diffuse fore- 
grounds rather than instrumental noise. In Fig. [T3] we 
show the comparison of the analyses which include the 
variation of the tensor-to-scalar ratio and the running 
of the spectral index, with fixed amplitudes of the point 
source contributions (red, dashed, curve) and without 
fixing the amplitudes of the point source contributions 
(black, solid, curves). 

Similarly to the standard ACDM model, leaving the 



Figure 12. Comparison of the cosmological and astrophysi- 
cal parameters obtained with the analysis which considers 
the standard De Zotti et al. based radio Poissonian (red 
curve) and the one which considers the Tucci et al. based 
radio Poissonian in the mock data and the standard De Zotti 
et al. based Poissonian in the code (black curve). 



astrophysical contributions free enlarges the errors bars 
including those of the running of the spectral index. The 
bi-dimensional plot of the scalar spectral index versus its 
running in Fig. ll4l shows. in addition to the enlargement 
of the errors, a degeneracy between the two cosmological 
parameters. The combined variation of the tensor-to- 
scalar ratio and running, without the inclusion of B- 
mode polarization shifts the mean of the running of the 
spectral index towards negative values with respect to 
input ones (n run — 0). 

We repeated also the main robustness analysis in 
this extended parameter space. In particular we re- 
peated the three main cases: the case with the Pois- 
sonian contribution multiplied by two, the variation of 
both amplitude and frequency dependence of the clus- 
tering term with the multiplication by two of the 100 
and 143 and multiplying and dividing by two the 217 
term. In Fig |15l we report the results of the analysis. 
We note how the recovery of the parameters is the same 
as in the standard astrophysical case showing how our 
method is robust also for this extended parameter space. 
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5 DISCUSSION AND CONCLUSIONS 

In the present study, we developed a multi-frequency 
approach of the analysis of the CMB power spectrum 
aiming at obtaining unbiased cosmological parameters 
in the presence of unresolved extragalactic point sources 
and TSZ signal. 

We have considered the three main extra-galactic 
contributions to the CMB signal at small scales, such 
as radio sources, infrared sources and SZ clusters. We 
have modelled these three signals in a minimal way us- 
ing the observational constraints from P lanck's early 
results (|The Planck Collaboration! l2011al lbh. In partic- 
ular, the radio and IR sources were modelled through a 
unique Poissonian term and a unique clustering term. In 
total, three amplitude terms suffice to account for and 
characterize the extra-galactic contributions. This min- 
imal number of parameters was obtained by considering 
the frequency dependences within our parametrizations 
of the TSZ signal, and IR and radio point sources. 

The use of a data-based physical model including 
the frequency dependence of the astrophysical signals 
under consideration to build the parametrizations of 
the small-scale contribution to the CMB is one of the 
original points of the prese n t work. Previous studies 
(|Efstathiou fe Grattonl (20091 ; iMillea et all l201lf ) chose 
blind parametrizations of the signals which consider 
generic shapes and parametrize each frequency with a 
different set of parameters. Such an approach applied 
for a frequency channels, leads to either a rapid in- 
crease of the number of paramete rs necessary to desc ribe 
the astrophysical signals, like in IMillea et al] (|201 J ), or 
it requ ires the use of a limited num ber of frequencies, 
like in lEfstathiou fc Grattonl HH). Other studies fo- 
cused only on one single contribution using more de- 
tailed and physical based parametrizations which in- 
volve a higher number of par ameters with respect to 
our ap proach (see for example lEfstathiou fc Migliaccirj 
(2011) for the TSZ) or considered the frequency depen- 
dence on a mu ch simpler leve l like for the clustering 
contribution in lAddison et al.l (|201lh . A different ap- 
proach, relying directly on the data, has been used b y 
the SPT team (|Keisler et al.ll201ll ; iLueker et aLlbOloT ). 
They took advantage of the high resolution to acquire 
data on scales where the CMB is suppressed. They used 
the data on smallest scales to measure the amplitudes of 
their blind templates for the astrophysical signals and 
consequently to remove the associated contamination 
from the CMB data. 

We have applied our multi-frequency minimal 
parametrization of the extra-galactic contamination to 
estimate the cosmological parameters for a standard six- 
parameter ACDM model and for a model with tensor 
perturbations and a running of the spectral index, with- 
out including residuals from our Galaxy and without 
using the B-mode information. We show that the best 
combination of frequency is the one which considers the 



frequencies from 70 to 353 GHz. It allows to estimate 
the cosmological parameters and at the same to put the 
tightest constraints on the amplitude of the IR cluster- 
ing term. We show, both in the standard ACDM model 
and in the "extended-parameter" model, that the input 
cosmological parameters together with the astrophysi- 
cal ones (amplitudes of the TSZ, Poisson and clustering 
terms) are recovered as unbiased and with slighlty en- 
larged error bars. We show that the Poisson and clus- 
tering terms are degenerate. 

Our approach instead of blind parametrizations 
prefers to use the whole available information on point 
sources and theoretical knowledge of the TSZ to have 
models of the astrophysical signals which are the most 
similar as possible to the expected signal. The use of 
the frequency dependence and the technique used to fit 
the data allowed a very restricted number of parameters 
to characterize the signals, in particular only a total of 
three parameters for all the signals. 

We have investigated for both cosmological mod- 
els, standard and extended, the dependence of our ap- 
proach and of the derived parameters on the accuracy 
of our model of extra-galactic contributions. In this way 
we tested its robustness against our theoretical priors 
on frequency dependence. We tested different Poisso- 
nian and clustering contributions by varying the asso- 
ciated terms both in amplitude and in frequency de- 
pendence. We showed that for both cosmological mod- 
els the cosmological parameters are recovered unbiased 
even in situations where the mock data differ by sev- 
eral sigmas from the model. The astrophysical parame- 
ters, amplitudes Aps and Acl, are recovered well in the 
case of a "real" Poissonian term different from that of 
the parametrisation. In the case of a different clustering 
term, the degeneracy with the Poissonian term induces 
a shift also in the Poissonian term. 

We have provided a minimal parametrisation to 
take into account the extragalactic contaminations on 
small scales. Our minimal analysis can be extended to 
include additional contributions like residuals from the 
Galaxy contamination. However, it is ideal in the sense 
that systematic effects, beam uncertainties etc are not 
accounted for. 
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Table 5. The exponents a, /3 and the coefficients A, B, C and 
w u for each term of the fits of the radio galaxies composite 
number counts. 
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Table 7. The coefficients and exponents of the fits for the 
clustering term Eq |2.9l 
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Table 6. The coefficients and exponents of the radio 
galaxy Poissonian contribution to the angular power spec- 
trum Eg [2741 



Figure 13. Results of the analysis which considers together 
with cosmological and astrophysical parameters also the vari- 
ation of the tensor to scalar ratio and the running of the 
spectral index. Comparison between the results varying or 
not the astrophysical contributions: the black (solid) curve 
considers the astrophysical signals whereas the red (dashed) 
curve represent the results with astrophysical contributions 
fixed to the fiducial value. Vertical bars are the input param- 
eters. 
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Figure 14. Results of the analysis which considers together 
with cosmological and astrophysical parameters also the vari- 
ation of the tensor to scalar ratio and the running of the 
spectral index. Comparison between the results varying or 
not the astrophysical contributions: the black (solid) curve 
considers the astrophysical signals whereas the red (dashed) 
curve represent the results with astrophysical contributions 
fixed to the fiducial value. We show the comparison of the two 
analysis for the running of the spectral index and the scalar 
spectral index. Curves are the 68% and 95% confidence level. 
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